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By Elbridge Z. St-well 


SUMMARY 


On the Basie of modern plasticity considerations, a unified 
theory of plastic Buckling appiicaBle to Both columns and plates 
has Been developed. For uniform compression, the theory shows 
that long columns which Bend without appreciable twisting require 
the tangent modulus and that long flanges which twist withoiit 
appreciable Bending require the secant modulus . Structures that 
Both Bend and twist when they Buckle require a modulus which is 
a combination of the secant modulus and the tangent modulus , 


INTRODUCTION 


The calculation of the critical compressive stress of columns 
and of structures made up of plates is an important problem in 
aircraft design. Formulas for the critical compressive stress 
have Been worked out for a multitude of cases of Both columns and 
plates. But these formulas are accurate only if the Buckling takes 
place within the elastic range of the material. In present-day 
designs, most Buckling occurs above the elastic range. The usual 
method of handling this problem is to retain all tire formulas 
derived for the elastic case. But to try to discover an effective, 
or reduced, modulus of elasticity which will give the correct result 
when inserted into these formulas. 

Column Buckling was the first structural problem to Be studied 
in the plastic range. In the latter part of the nineteenth century, 
Engesser proposed use of the tangent modulus as the reduced modulus 
for columns. At almost the same time, in the Belief that- the 

column would Be strengthened By unloading on the convex side, 

Considers suggested that the effective modulus should lie Between 
the tangent modulus and Young* s modulus. This concept was subse- 
quently refined By Engesser and By von Karman (reference l) and 
led to what is generally known as the "double modulus . " 
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Experiments have shown, however, that the von Karman double 
modulus gives values that are too high tor the column strength 
(reference 2) and that the correct modulus is nr oh ably the 
tangent modulus. Shanley (reference 3) has stated the situation 
compactly as follows: "If the tangent modulus is used directly 

in the Euler formula, the resulting critical load is Bomewhat 
lower than that given by the reduced modulus theory. This 
simpler formula, originally proposed by Engesser, is now widely 
used by engineers, since It gives values that agree very well 
with test data." Further careful tests by Shanley (reference 4) 
and also by Langley structures research laboratory have shown 
that the unloading on one side of the column, postulated by 
von Karmen, does not occur at buckling rnd that the correct modulus 
for columns is actually the tangent modi’. !.us. This conclusion also 
has theoretical justification (references 3 and 4). 

In. the case of local or plate buckling, the reduced modulus 
is appreciably higher than the tangent modulus. Tests of the 
local budding stress of aircraft— section columns have been made 
by Gerard (reference 5), who has suggested the use of the secant 
modulus for this type of buckling. Extensive tests in the 
Langley structures research laboratory on similar aircraft sec- 
tions made and reported over a period of several years and 
summarized in reference 6 have also shown that the reduced 
modulus for plates is in the vicinity of the secant modulus. 

In particular, tests of long alumj.nuEv-a.lloy cruciform-section 
columns, designed to buckle by twisting without appreciable 
bending, have been made in a manner similar to that described 
for the aircraft— section columns in reference 6. The results 
have shown that the reduced modulus for pure twisting is very 
close to the secant modulus. 

The present paper constitutes a theoretical investigation 
of the buckling of plates beyond the elastic range, which includes 
columns as a limiting case. Such an investigation requires a 
knowledge of the relations between, the stress and strain components 
beyond the elastic range . These relations have not as yet been 
conclusively determined. A recent paper by Handelxnan and Proger 
(reference 7) based on one possible set of stress— strain relations 
led to results for -the buckling -of hingod flanges in sharp disagree- 
ment with test results obtained at the Langley structures research 
laboratory. Another set of stress-Mstraiu relations is generally 
accepted by the- Russian investigators and has been, applied by 
Ilyushin (reference 8) to the stress conditions in thin plates. 

These results form the foundation of the present paper, which 
assumes that in plates as Well as in columns, unloading during the, 
early stages of ■ buckling does not occur. On this basis, a unified' 
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theory of plastic "buckling applicable to both column and local 
buckling has been developed. The results are presented in the 
following section. 


RESULTS ML CONCLUSIONS 


Ilyushin (reference 8 ) has treated the stability of plates ■ 
stressed above the -elastic limit with consideration of the three 
possible zones that might result from buckling: ( 1 ) a purely 

elastic zone, ( 2 ) a zone in which part of- the material is in the 
elastic and part is in the plastic state - the "eXasto-plastic" 
zone, and ( 3 ) a purely plastic zone in which all of the plate 
is stressed beyond the elastic limit. All three zones may exist 
simultaneously if the plate is not entirely in the plastic state 
before buckling or if the buckling is allowed to proceed beyond 
the initial, stages. 

If, however, the plate is uniformly loaded before buckling 
so that all parts of it aro initially at the same point in the 
plastic range and if, in addition, buckling and increase in load 
are assumed to progress simultaneously, then the plate may be 
expected to remain in the purely plastic state in the early stages 
of buckling. This second assumption is in agreement with the 
corresponding condition that apparently holds for columns 
(reference 4). 

Upon the assumption that the plate remains in the purely 
plastic state during buckling, Ilyushin's general relations for 
this state have been used to derive the differential equation of 
equilibrium of the plate under combined loads. Since critical 
stresses are generally simpler to compute from energy expressions 
than from a differential equation, the corresponding energy expressions 
were also found. These derivations are given in appendix A, 
together with applications to compressive buckling of various . types 
of plates. A comparison with Ilyushin's treatment of the plastic- 
buckling problem is given in appendix B. 

The results of most interest in the present analysis are given 
in the following table as values of a quantity q, the number by 
which the critical stress computed for the elastic case must bo multi- 
plied to give the critical stress for the plastic case: 
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These values of r\ are plotted as curves A to G in figure 1 
for extruded 24S— T aluminum alloy for which the compressive yield 
stress was 46 ksi. Similar curves for rj could readily "be 
prepared for any other material having a known stress-strain 
relationship. 
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The values of t) given in the table were obtained by dividing 
the critical stress of the structure in the plastic region by the 
critical stress that would be obtained on the assumption of perfoct 
elasticity. Since Poisson's ratio has been taken as one-half in both 
computations f errors ffom this cause will ordinarily be present in 
both critical stresses. Most of these errors will be eliminated, 
however, in the process of division to obtain and, consequently, 
the values of rj given are believed to be nearly correct. 

When plate -buckling stresses in the plastic range are to be com- 
puted, the experimental value of Poisson's ratio' that applies as 
closely as possible to the stressed mate ial, together with the appro- 
priate value of T) from this paper, should be used in the plate- 
buckling formula. 

The highest value of rj which is . can be realized only if 

E 

there is negligible longitudinal bending (as with a long hinged flange 

E + 

which buckles by twisting) . The lowest value of rj which is van 

‘ •• E 

occurs when the longitudinal bending predominates over other types of 
distortion (as with a long column under Euler buckling). The theory 
implies that a change in the stress -strain curve caused by pre- 
stressing of the material, would alter the value of tj in the first 
case but not in the second, if the buckling stress is higher than the 
highest stress reached during the opera tiqn of pre-stressing. If, cn 
the other hand, the buckling stress is lower than the highest stress 
reached during the operation of pre -stressing, then = 1 for each 
case. 


National Advisory Committee for Aeronautics 

Langley Memorial Aeronautical Laboratory 
Langley Field, Va., July 2$, 19^7 
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APPENDIX A • ■ 

THEORETICAL DERIVATIONS 

Dof jnltions .— The intensities of stress and strain are defined 
in reference 8, respectively, as 



where 

stress in the z-dir action 

€ x strain in the x-direotion 

Oy. stress in the y-direction 

strain in the y-direction 

t shear stress 

7 sheen 1 strain 

According to the fundamental hypothesis of the theory of plasticity, 
the intensify of stress is a uniquely defined, single— valued 

function of the intensity of strain e^ for any given material 

if increases in magnitude (loading condition). If 

deoroases (unloading oondition), the relation hetwoen a ^ and e^^ 

hecomes linear as in a purely elastio case. 

In the equations of definition (l) and (2), the material is tahon 

to ho incompressible and Poisson’s ratio = i. The sti'oss— strain 

2 

relations compatible with the equations of definition (l) and (2) are: 
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6 

X 



E 


sec 



E 

eec 


e = 
7 


°7 ' t* 1 s y 


E, 


sec 


E 


sec 


> 


7 » 


3t 


E 


sec 


e i = 




E 


'sec 


(3) 


These relations imply isotropy of the material . 

Variations of strain and stress * - When buckling occurs, let 
G yt 311(1 7 vary slightly from, their values before buckling. 

The variations Be x , 5 e y , and 67 will arise partly from the 

variations of middle- surface strains and partly from strains due to 
bending j thus. 



Z *1 



> 


(*) 


57 = 26 ^ - 2zx^ 


in which and e 2 are middle -surface strain variations and 

is the middle-surface shear-strain variation, and Xg are the 
changes in curvature and X^ is the change in twist, and 
distance out from the middle surface of the plate. 


z is the 
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The corresponding variations 5S X , and 5t in 

S v , and t must he computed. From equations ( 3 ) > 


S_. “ € 

x sec x 


therefore 


5S X ^ E seo 5{? x + e x 


c 6 f~ N ) 

w 


= E sec ^ e x 


_ 1* _ 

e i \ «i 


do, . 

dej 1 ■*• 


(5) 


Wow the variation of the work of the internal forces is 


a l “ °x ^- € x + °y ® f y + 


so that 


5ej 


a x 8e x + a y 6e y + T 


a x € l 


y 2 


(.Vi 


+ or + 2 Te„ - z(o x, + a + 2xX 




Si 


( 6 ) 


Substitution of this value of 5e^ in equation ( 5 ) gives 


5S x = E sec 56 x ~ 


i*. Al _ fV 

ff i°i V 1 d0 i> 


a X e l + V 2 + ^*3 


“ z ( a x\ + a y\ + 2tX 3) 
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Let the coordinate of the surface for which 8e^ = 0 (the neutral 
surface) he z = z Q . The expression for z Q is obtained by setting 

8e^ =0 in equation (6); 


z 


o 


p x g i + Va + 
tT x Xi + cr y x 2 + 2 t^ 3 


(T) 


By introduction of this coordinate into the expression -for 8S X and 

p i dcf^ tu 

by recogjiition of — as EggQ deT aS tan 


SS = E 
x sec 


(fl ~ zX l) 


+ ^ E sec ~ E tan. 


ffiSiV 


) (Vi + p y X 2 + 2tX 3 ) < z “ z ° 


) ( 8 ) 


In a similar way it may be shown that 


5S y = E sec 


°i e i 


(*2 - « 2 ) 

( E eao - E tan) (vi + Va + 2rX ^ Z ~ Z ° ] (9) 


and 


8r = — E 


see 


3 a i e i 


(^3 " ZX 3 ) 

( E sec " ^tan) ( p x*l + Va + 2tX 3 ) < z “ z o> < 10 ) 
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Variations 
Impressed forces 

and 


’ forces and moments .— For the variations of the 

T y) T , and T_,,_ and the moments M.., M . 

A o -Of «y 


5 T_ 


ST 


f | 

1 / l Scr x dz 
*• 2 - ■ •• 

h 

h 5c V dz 

0 2 

r.& 


8T 


xy 


u 


h 

‘2 ' 

h 


5t dz 


> 


6^ = = 5 d x z dz 

U‘2 ’ • 

ft 


SM 


y k y 

'-'“2 


Scr z dz 


6M * 


it 

“S 


8t z dz 


( 11 ) 


where h is the thickness of the plate 
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where 
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Similarly, 



( 14 ) 


In these expressions,, the Integrations of 6S , 6S , and 8t in 

. x y 

the plastic region have heen taken over the entire thickness of the 
plate, with the Implication that no part of the plate is being 
unloaded. ... . . 


Bcmatlon of equilibriu m.— If w(x,y) is the bending deflection 
of the plate at buckling, and if no external moments are applied to 
the plate, then the equation of equilibrium of an element of the 
plate may be written 


bx 2 


+ 2 


a2 ( 6M y) 


bx by 


by* 


('* 


q 2 w 


dx' 


2 + ^S +27 


cfir 
CSX by. 


( 15 ) 


in which the impressed forces c^h, a h, and ih are considered 
as given (cr x and cr y are positive for compression). In terms of 
w, the changes in curvatures are 
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and 


The change in twist is 


X 

2 


b%r 

Sy 2 


X 

3 


5 2 W 
cbc dy 


(3.6a) 


(i6h) 


( 16 c) 


When the values of SM^., SMy, and SI- r y in equations (12), (13), 

and (14) , respectively, are differentiated as required by equation ( 15 ) 
and substituted in that equation, the general differential equation 
of equilibrium for a plate in the plastic state is obtained as follows; 



(17) 


In the elastic range, equation ( 17 ) reduces to the usual form 

■s 
02 


vV = - | 




+ 2t 

Sy" Sr oy 



where 
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Energy integrals.- Equation (17) is the Euler equation that 
results from a minimization of the integral 





Six 2 bx by 


, /^w\ 2 

3 s 2 * Sjr + c / ^i 

^ Sx by c5y 2 ^ V^y 2 / j 2 


/ \ 2 t b 2 * S 2 ^ 

\5x 5yy Sx 2 By 2 


fx /&w \ 2 

ff i vW 


+ £l&|e + 

ffj_ ox oy cr 1 \byj 


dx dy 


(18) 


/ 


which represents the difference between the strain energy in the 
plate and the worlc. done on the plate by the external forces. The 
coefficients in this ‘integral are : 


In the plastic region 

In the elastic region 

_ / <y \ 2 / ' E \ 

c, a 1 - l(~ ) (1 _ _£ss. } 
1 4 Vo, j \ bj } 

N 1 \ sec/ 

0 !-l 

CT 1 \ ^'sec/ 

c 2 = 0 

■n. -i_3 V.v + 2t ® ( E tan \ 

C 3 « 1 

’ 

3 4. \ ~) 

1 \ see/ 

cr T / I. \ 

0 4 o 3 -% 

i \ . sec/ 

0 

n 

-=* 

0 

, /cr_\ 2 / E \ 
5 4 W V E sec/ 

c_ =» 1 
5 


€ 









If there is a restraint of magnitude e along one longitudinal edge of the plate, the strain energy 
in this restraint Itself is taken to he 



(19) 


if j Q is the edge coordinate. (See reference 9 for form of expression,) In expression (l9),the 

stiffness D 1 is assumed to he the same as that in equation (12) . If restraints are present along 
two edges, there will he two terms similar to expression (19) . These terms may he added to 
integral (18) aB additional strain energy. 

Critical stress in ■nlaatic region .— If the integral (l8), supplemented if necessary hy addi- 
tional terms of the form of expression (19), iB set equal to zero and tbs resulting equation solved 



. . , . ■ -i* ‘ (•* ’ ' 

in which the values of the C’s in the plastic range are used. This expression for the critical- 
stress intensity may he minimized as with the corresponding elastic case , 


H 

VJl 
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If the values of the C*s In the elastic region are used In formula (20), the 
critical-stress intensity in the elastic region (di) e q la as follows: 


W) 


D 

el " h 


'i£\ , / aPr Y, && , (m 

sSk I. 2 / ■ Is® dyy 5x By 2 V^j 2 / 


dxdy + — 

D 


/Sw 

W 




dx 


y=y 


.al 


— (^f + 2 — Ididy 

a i VSx / cr^ 5x 3y , 


^xpresclon for n .~ A quantity q is defined as 

W, 


n 53 


c°i) 




el 


(21) 


(22) 


This quantity is a direct measure of the effectiveness- of p' nsticity in reducing the 
critical stress of a structure, and its computation in terms of the constants of the 
stress-strain curve represents the solution of the problem of plastic buckling. 


Application to platen compressed in the x-direction .— The theory will now be' 
applied to flat, rectangular plates uniformly compressed in the x— direction. Values 
of q will be computed for the following cases: 


I. Long plates with one free edge (flanges), the other edge being either hinged 
or clamped 


H • Long plates with both edges either hinged or clamped 
HI- Plates with two free edges (columns) 


o\ 
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When (T =s t w Oj a- » a and the plasticity coefficients reduce to 

J . * A 




c 3 L 3 _ta£t 
L 1 S * 4 E . 


QQC 


and 


Cg = =5 0 


°3 " °5 ' 1 


The differential equation of equilibrium, equation (17) , than "becomes 


„ aV : „ aV . a 4 » H 

1 a** a* 2 a^ aj 1 * »' 


cftr 


(23) 


and' the corresponding' energy' expression (20) for tlie critical stress in the plastic range 



Case I: Elastically restrained flange 


If y * 0 is the elastically restrained edge of the flange and y *» "b is the free 
edge, a deflection surface knovn to he good in the elastic range and presumably satis- 
factory also "beyond this range is (reference 9) 
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w = 


H 

00 


where 


r 




= -J+.963 

a 2 = 9.852 
a 3 = - 9.778 


and e is the magnitude of the elastic restraint. Substitution of this expression for v 

2 


^■7 


L_> — ’ — “ 

L _ I o,\ + i! 

1 c 5 

A + v + 

V s 2 y k 

c 6 2° 7 + irbN 2 

L %v . 

^ ^ 8 J'i) 


«) 2 A« 


1 c 8 € 
t + — + 
3 2a^ 


where 


c n e2 


4 a 2 

3 


b 2 h 


C 1 = 0.23694 

C 2 = 0.79546 
c 3 = 0.89395 

= 0,04286 


c 5 = 0.56712 
Cg = 0.17564 
c ? ? 0.19736 
Cg - -2.3168 

4.0982 


pi 
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In order to find the minimum Tain© of 



0 


vhich give a 



The minimum value of 





is therefore 



For the elastic case, the same expression is obtained from agnation (21) with = 1 and 
D T replaced by D. From equation (2£), therefore. 



H 

VO 
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This valy©., as ^function of stress, is plotted as curve A- for aluminum alloy in ; 
figure 1. The individual points represent the NACA tests of the buckling of cruciform- 
section columns for which' the condition € - 0 is fulfilled. 

(b) If the edge y = 0 1 b clamped, € = «° and, from equation (25), 
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This value of T) is plotted as curve B for 24S-*E‘ aluminum alloy 
in figure 1. 

Case IX: Plate elastically restrained along two Unloaded edges 

If y s + | are the immovable unloaded edges which are 

elastically restrained against rotation by restraints of magnitude 
satisfactory deflection surface is known to be (reference 10) 



Substitution of this expression for w in equation (21}-) gives 



where 

0.0237S 2 + 0.297 e + ~ 

.*!(«>- 

0.0046le 2 + 0.09k7e + - 
, 2 
and 

, , 0.01l4e 2 + 0.l89ke + 1 

fg(e) = * • — 

0, 00461 e 2 4- 0.0947€ + i 

2 

In order to, find the minimum value of ( a \ 



which gives 
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The minimum value of 



is therefore 



2 \fi * r a (c) 


rt 2 D r 

b 2 h 


For the elastic case, the same expression is obtained from 
equation (21) with «= 1 and D* replaced by D. From 
equation (22), therefore, 


E 2\j£J7) + f a (0 


(28) 


(a) If the edges y = t ~ axe hinged, e ■=» 0, f^(e) a 1 , 

f 2 (e) = 2, and, from equation (28), 

E___ 1 + \/c~ E 


7 ] = 


' sec 1 + V C l 
E o . 


E. 




sec 


1.1 1 + 3 . 

~ T — \ TT * \ 11^ 


IT \2 £\ F 4 E 


(29) 


sec 


This value of q is plotted as curve C for 243— T aluminum a l l oy in 
figure 1 . 


(b) If the edges y *» t -- are clamped, € = °°, f-j^e) “ 5*15* 
fp( e) = 2.46, and, from equation (28), 


E 


eec 

E 


2.46 + 4.52 ,/c i 


6.98 


E ' 
sec 

E 


,.353 + 0.64^| + f^2) (30) 


This value of q Is plotted as curve D for 24S— T aluminum alloy 
in figure 1. 
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Case III: Plastic buckling of columns . ' • 

For the discussion of the plastic buckling of columns, it is 
convenient to revert to the differential equation (23). The plate, 
when loaded as a column, has two. free edges described by the condi- 
tions 


£3 + I £z\ =0 

by 2 2 by 2 Jy=^& 



3 \ 

2 ox 2 SyJ+b 
~o 


A solution of equation (23) which identically satisfies the .first' 
condition is 

/ • ’ V . ... \ 

w = ( q cos — cosh — + p cosh — cos ^ j cos — 

2 b 2 b / l 


where 






(?*)pl " k 


b 2 h 


1 ■ + 5 if/ 


sec 
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In order that this solution also satisfy the second condition 
at the free edges, it' is required that 


a 2 q 


-0 


tanh ~ 

2_ 

a 

2 




Q. + 


(?) 


tan 

“T 

2 


(3D 


which is the buckling criterion for the plate when loaded as a 
column, let 


k + (i) a ( 1 -°i)’(i) a ( 1 -‘ a ) <3£) 

p 

where 6 Is a quantity to be determined for three individual cases. 
By use of equation (32) ' ' 

a 2 = (l + \[l - £ 2 ^) 

"* ■ (r)^ 1 + N ^ 2 ) 

’ - (f T (l * vRj) 

s = (f ) 2 (- i * x / 1 - 5£ ). 

/nbf 

P- 5 =( r ) 

and the buckling criterion given 


in equation ( 31 ) becomes 


2 2 tan I 

f3 2 p 2 ± = 0 

* £ 
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or 




1 


tcnh | 


a, 

2 


i 2 + 


tan — 


£ 

J 2 


= 0 ( 33 ) 


y'% 

From equation (32), k = ^ ^ - £ 2 ^; and thus the critical 

(<w3. 


stress in the plastic range is 

(%) 


jt 2 E 


sec 


'pl 


i © 5 


where 


P = 


\/l2 


The corresponding critical stress in the elastic range is 

, ^ * 2 E(l - f~) 

w “ ' ~W~ 


The reduction factor q is obtained from formula (22) as 


(M 


E 


sec 

E 


E 

3 tan 
¥ E 


( 34 ) 


1-1 


(a) In order to investigate the case of short columns, let £ 
approach” zero. Then, by definition of 3 ,' ’ 
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. and 


tan- 

T 


£ 


2 


In addition, 

i 2 * (| - S 2 ) (l - ^/l - 5 2 )— >0 
The buckling criterion given in. equation ( 33 ) therefore 



1 

The expression in the brackets approaches ~ as §■— >0 
to satisfy the buckling criterion, therefore. 


tanh — 


— ; 

£ — >0 

a 


2 

which can be realized only if 
For short columns, therefore, 

a is large j that is, if 


S » 0 

and, from equation (34), 

• 


E ' E 

1 - f + a -ifs 

4 e 4 E 


reduces to 


In order 


— is large. 

I 


(35) 


This value of q is plotted as curve E for 24S— T aluminum alloy 
in figure 1 . 
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("b) For a squar e plate , -. = 1, a = it\jl + \Jl — i 2 , 

P = lie \J± - sfl — and the "buckling criterion given in 

equation (33) "becomes 


S 2 + (I - ?) (*■ + i 1 - 5 s ) 


tash(j\/l <-^1 - 8^ 
|\/l +N /l- 6 2 




tanh 

V 1 

'1 -i s ) 

sJ 


-£ 2 


which is satisfied "by 5 = 0,15375. From equation (3k), 

n - 0,11k ?222. + 0.886 

1 E 


(36) 


This value of q is plotted as curve F for 2kS— T aluminum alloy 
in figure 1. 


(o) For long columns, a and £3 "become so small that 

a 
2 


tanh — tan — 


a ■ . J3 
2 -2 


4«i 


and the "buckling criterion of equation (33) reduces to 

s - *§~ s j i 1 + n / 1 - ^ - [ ?2 + f) -f- 


p n 

which is satisfied "by | s i, From, equation (3k), 


E, 


tan 

33 


(37) 


This value of rj agrees with the experimental results of 
references 2 and k and is plotted as curve G for 2kS~T aluminum 
alloy in figure 1. 
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APPMUEX B 

COMPAEESON HITH ILYUSHIN'S STABILITY CALCULATIONS 


The basic difference between Ilyushin's solution of the plastic- 
buckling problem and that given in this paper is that Ilyushin con- 
siders the plate to unload on one face as it-buckles. The unloading 
process reaults in the creation of an elasto-plastic zone in the 
plate, and different equations from those that apply when the plate 
remains plastic during the buckling process are required for this 
zone. 


The differential equation for the buckling of a rectangular 
plate when buckling is accompanied by unloading is given by Ilyushin 
as equation (3.43) of reference 8. For simple compression in the 
x-direction this equation is of the same form as equation (23) of 
the present paper, but with the following different constants: 

D is used instead of D' and 

k =» 1 - l£ 2 (3 - 2$) 

E 

is used instead of — , Tn the formula for k, from equation (3#l) 

®sec , 

of reference 8, 


i 



and, from equation (It 22) of reference 8, 



E 


■When the values of £ and \ are inserted into the expression 
for k 
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E 

Computation shows that k Is always larger then , so that the 

E tan EseC 

use of k in place of will result in appreciably higher values 

®sec 

of r] than those given in the present paper. Since IlyuBhin uses 
the elastic value D * there is no possibility of the solution ■ 
yielding a secant modulus. Curves ' A to 0 in figure 1, if com- 
puted from Ilyushin’s equation (3»^3)> would start with a horizontal 
line at unity for curve A (Young's modulus) and end with curve G 
expressing the Karman double modulus which is appreciably higher than 
the tangent modulus of this paper. If D’ were substituted, for D 
in Ilyushin's equation (3«^3 ) } curve 1 A would then represent the 
secant modulus a3 it does in the present paper, but curve G would 
still remain the Karman double modulus. Therefore, when the unloading 
of the plate during the buckling process is considered, results are 
obtained which are not confirmed by experiment. 
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Figure I. — Computed curves showing variation of r/ with stress for 
various structures of 24S-T aluminum alloy in compression. (Curves 
A to G are drawn for a material with a yield stress of 46 ksi. The 
scatter band enclosing curve A shows the limits of variation of specimen 
properties from this value.) 



